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Abstract 

We investigate a propagation of solitons for nonlinear Schrodinger 
equation under small driving force. The driving force passes the res- 
onance. The process of scattering on the resonance leads to changing 
of number of solitons. After the resonance the number of solitons 
depends on the amplitude of the driving force. 

Nonlinear Schrodinger equation (NLSE) is a mathematical model for wide 
class of wave phenomenons from signal propagation into optical fibre PI E] to 
surface wave propagation [Hj . This equation is integrable by inverse scattering 
transform method |1] and can be considered as an ideal model equation. The 
perturbations of this ideal model lead to nonintegrable equations. Here we 
consider such nonintegrable example which is NLSE perturbed driving force. 

The most known class of the solutions of NLSE is solitons The struc- 
ture of this kind of solutions is not changed in a case of nonperturbed NLSE. 
The perturbations usually lead to modulation of parameters of solitons IH] . 
Number of solitons does not change. 

In this work we investigate a new effect called scattering of solitons on 
resonance. We consider the process of scattering in detail and obtain the 
connection formula between pre-resonance and post-resonance solutions. In 
general case the passage through resonance leads to changing of the number 
of solitons. This effect is based on the soliton generation due to passage 
through resonance by external driving force |7|. 

We found that the scattering of solitary waves on resonance is a general 
effect for nonlinear equations described the wave propagation. In this work 
we investigate this effect for the simplest model. It allows to show the essence 
of this effect without unnecessary details. 

iThis work was supported by RFBR 03-01-00716, grant for Sci. Schools 1446.2003.1 
and INTAS 03-51-4286 
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This paper has the following structure. The first section contains the 
statement of the problem and the main result. The second section contains 
the asymptotic construction in the pre-resonance domain. In the third section 
we construct the asymptotic solution in the neighborhood of the resonance 
curve. The fourth section of the paper is devoted to construction of the 
post-resonance asymptotics. Asymptotics are constructed by multiple scale 
method jS] and matched jU]. 

1 Statement of the problem and result 

Let us consider the perturbed NLSE 

idt^ + dl^ + 1^1'^'^ = e^fe'^^'\ < e < 1. (1) 

The phase of the driving force is S/e"^ = ut. The amplitude / = f{ex) is a 
smooth and rapidly vanished function. 

In the simplest case the phase is linear function with respect to t S/e'^ = 
ut, uj = const. In general situation the constant frequency of the driving 
force does not lead to scattering of solitons. Let us investigate the driving 
force with slowly varying frequency. The most simplest dependence on t for 
UJ has a form u = eH/2. The amplitude / of the driving force admit an 
additional dependency on eH but it leads to complicated formulas and no 
more. 

Let us formulate the result of this work. Below we use the following 
variables Xj = e^x, tj = eH, j = 1, 2. 
Let the asymptotic solution of ((H) be 

"^{x^t^e) = e u {xi,t2) + Oie"^) as ^2 < 0, 

where u {xi,t2) satisfies 

dt2 u +d^^ u+\u\^ u=0 

and initial condition 

u \t2=to = hi{xi), t° = const < 0. 
Then in the domain ^2 > the asymptotic solution of has a form 



"^{x, t,e) = ev (xi, t2) + 0{e^), 



(2) 
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where v {xi,t2) is a solution of NLSE with initial condition 

V |t2=o =u {xi, 0) + (1 - i)^/^rf{xl). 

Let us explain the result for soliton solution. If in the domain ^2 < the 
solution has A^-soliton form then in the domain ^2 > the number of solitons 
is defined by initial condition for v. 



2 Incident wave 

In this section we construct the asymptotic solution of equation in pre- 
resonance domain. This solution contains two parts. The first part is a 
specific solution of the nonhomogeneous equation. This solution oscillates 
with the frequency of the driving force. The amplitude is determined by 
an algebraic equation. The second part of the solution is a solution of the 
homogeneous equation. The solution contains undefined function due to 
integration. This undefined function usually determines by initial condition 
for Cauchy problem. 

We construct the formal asymptotic solution of the form 

132 
-^{x, t,e)=eu (xi, t2) u (xi, ^2) + B (xi, ^2) exp(iS'/£^) + 

e^l^Bi {xiM)^w{iS/e^)+ B-i (xi, t2) exp(-i5/e2)^ (3) 

B2 (xi,t2)exp(2i5/e2). 

To determine the coefficients of the asymptotics substitute (jH} into equa- 
tion (H)). It yields 

e^i^ - S' B -/^ exp(iS'/£:^) + e^[i + u^^^^ +| m |^ u ^ 
+^^{f^- S' Bi +i Bt2 + Bx^x^ +2| B^ exp{iS/£^) + 

+ (^S' B-i + h^B*^ exp(-z5/£2)^ 

3 3 .,ili2 3 I23* r,l|2.o 

e \lUt^+ Majixi +2| u \ u + u u + 2u \ B \ + 



-2S' B2 + U* B^\ exv{2iS/e^)] = 0(£^i?(t2, xi)) 
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The residue part of the asymptotics has a form 

R{t2,x,) = 0{\ B + u \^ + e^\ Bi I' + e'l B-i I' + e'l B2 ?)■ 

Collect the terms with the same order of e up to the order of and 
reduce similar terms. It yields differential equations for m, u and algebraic 

2 4 5 

equations for B, B±i and B2- 

.1 1 I 1 12 1 r^ 

I Ut^ + U^^^^ +1 M I U=0, 

• 3 3 „i 1 i2 3 1 2 3 * r^i 2 ,2 1 

t + Mzixi +-^1 U \ U + U U = —2\ B I U, 



- S' B= /, 

-5' Bi= I Bp, + Ex^ixi +1 u |2 5, 
-25' 52=^^ * B 



B-i=-u^B\ 



1 3 

The coefficients u, u are uniquely determined by initial conditions at the 
moment ^2 = ■ We suppose that = const < and 

1 3 

where functions hi , are smooth and rapidly vanish I — s> ±00. 

The coefficients of the representation Q have a singularity as S" — > 0. 

The order of singularity of Bk is easy calculated. 

B= 0{t-^), Bi= 0{t-'). 
3 

To determine the asymptotics of m as ^2 ~0 we construct the solution 
of the form 

u= t^^ u (-i'°)(xi,t2) + ln|t2| u ^^'^\xi,t2)+t2\n\t2\ u ^^'^^ (xi , ta) + ^(a;! , ^2) ■ 



3 

Substitute this representation into equation for u and collect the terms of 
the same order with respect to ^2- It yields equations for coefficients of the 
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asymptotics 

3(-i,o)^^2|/|2i, 
u (O'l) = -iL{u (-1-0)), 

L{u) = it2\n\t2\L{l ^'^'y) + I u (4) 
Here ^^(m) is a linear operator of the form 

L{u) = idt^u + dl^u + 2\u pu+ m 

Functions U (-I'O), M (0,1) 

and u (^'^) are determined from algebraic equations. 
These functions are bounded as —const < ^2 < 0, const > 0. 

3 

The function m is a solution of nonhomogeneous linearized Schrodinger 
equation. The right hand side of the equation is a smooth function as 
—const < ^2 < 0, const > 0. The solution of this equation can be ob- 
tained using results of jTO] . In particularly if u is N-solitons solution of NLSE 
then exists the bounded solution of nonhomogeneous linearized Schrodinger 
equation (jU as —const < ^2 < 0, const > 0. 

Coefficients of Q have singularity at t2 = 0. After substitution Q into 
equation we obtain a residue part. This residue part increases as t2 ^ 0. 
The domain of validity of (jSI) is determined by 

e^R{t2,xi) = o{e), e ^ 0. 

It yields 

-t2>£: or -t>£:"\ 

3 Scattering 

In the neighborhood of the point t2 = the frequency of the driving force 
becomes resonant. It leads to changing of behaviour of the system described 
by (0). The forced mode of oscillations changes by a resonant mode. Formally 
it means representation Q is not valid. 

In this part of the work we construct another representation for the solu- 
tion of equation (Q). This representation is valid in the neighborhood of the 
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resonance line ^2 = 0. 

1 2 

\l/(a;, t,e) = e w {xi,ti)) + w (xi, ti) + 

3 1 3 

e^lne w (xi,ti) + w {xi,ti) e ^ 0. (5) 

Here we use a new scaled variable ti = e. Representation (0) is matched 
with (jni)- It means these formulas are equivalent up to value o{e^) as ^2 — 0. 
The coefficients S of © are determined by ordinary differential equations 
©, (IHl), (fTIH) and matching conditions. 

To obtain the behaviour of the coefficients of (0) as ti — > — cxd match (0) 
with (jS)). Write (jSJ) in terms of ti 

=el^u (xi,0) - {t^^ f + it^^ f) exp{iS / e^)^ + 

e^{dt, u {xi,t2)\t,=oti+t^'i\f\^ u (xi,0) +0(tr')) + 
eHne^-iL{2i\f\^u)\t,=, + o{l)^ + 

+^^(^^'2 ^ (^i'^2)|t2=ot? + ^(a;i,0) + o(l)^, e^O. 

To obtain equations for coefficients of ^ substitute (jS)) into equation ((T)). 
It yields 

e^[{dt^ w -f exp{iS/e^)^ + e^(^dt^ w +dl^ w +7I w + 
e'^(^dt^ w +dl^ w ++ w'^ w* + 2-^\w = 0{e^p{ti,xi,e)). 

The function p(ti,Xi,e) can be represented in the form 

xi, e) = 0(1 w p w +dl^ w +e\ w \^ + e'^l w f). 

Collect the terms of the same order with respect to e. As result we obtain 
the equation for w 

tdt,w= fexp{itl/2). (6) 
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The matching conditions give w=u (a;i,0) ti — oo. The solution of this 
problem is represented in terms of Fresnel integral 



w=u {xi,0) - if (xi) r exp{iey2)de. (7) 

J —oo 

Equations for higher-order terms are 

idt-^ w= —dl^ w +\ w l"^ w, (8) 

tdt, w= 0, (9) 

d3 2 |l|o2 1 o 2 ^ / -1 r\\ 

w= —o^^ w — 2| w \ w — w w . (10) 

The higher-order terms satisfy fist order ordinary differential equations 
with respect to ti. The spatial variable xi is a parameter in these equations. 
The solutions of these equation are uniquely defined by terms of the order of 
1 in asymptotics as ti —>■ —oo. The asymptotics as ti — > — cxo is obtained by 
matching 

2 1 

w= dt^ u (xijta) 1*2=0^1 + 
^i = -^L(2z|/pi)|,,=o + o(l), 



To determine the behaviour of the solution after resonance we need to cal- 
culate the asymptotics as r — *• -|-oo of the coefficients for representation (jH)). 
Calculations give 



w (xi, ti) =u (xi, 0) - if{xi) 



.^^ _^ exp(a;/2) _^ 



where Ci = (1 — i)\fn. 
Denote by 



2 

The function w has the asymptotics of the form 

w (xuh) =tiw & o(xi) + + O(tr^) 
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where 



^ 1 = -^xi ^0 +1 w^o I 




dii^i) = kidx^f + ^1 and are constants. 



w {xiM) =t\w 2{xi) + o{t\). 



where 



w 2(3:1) = I [d^^ wi +2\ wq y wi + w Q w ^ 



Representation (0) of the solution for (P) is vahd as 

e^p{ti,xi,e) = o{e). 

The determined above behaviour of coefficients of asymptotics ^ give the 
domain of vahdity with respect to ti 

< E^'^ or \t\ < 

4 Scattered wave 

In this section we construct the asymptotic solution of equation after 
the resonance. The leading-order term of the solution satisfies NLSE and 
depends on Xi,t2 as well as before resonance. But this leading-order term is 
determined by another solution of NLSE which contains generally speaking 
another number of solitons. This number depends on a condition on the 
resonance curve ^2 = 0. 

After resonance we construct the asymptotic solution of the form 



*(a;, t,e) = e V {xi,t2) + e'^ v (xi, ^2) + 
A (t2, xi) exp{iS/e^) + e^{Ai (^2, xi) exp{iS/e^) + 

A~i (t2,xi)exp{-iS/e )). 



(11) 



Substitute this representation into ((H): 



e^-S' A -/) expitS/s^) + e^dt, I +dl 



1 I 1 i2 1\ 
V +\V \ V) + 
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e'^{dt^ V +dl^ V +2\v l"^ V + v V* + 
i-S' Ai +dt, A +dl A +2| h \^ A) expiiS/e^) + 
{S' A-i + h^A*) exp{-tS/e^)) = eh{t2, xi,e). 
Here r(t2, 3^1,^) depends on coefficients of the asymptotics (fTT|) . This depen- 

2 4 4 

dence is easy calculated. The coefficients A, Ai A-i have singularity on 
the resonance curve. To determine the domain of validity of we need to 
derive the explicit formula for r 

r{t2,xi,e) =0{l + e\ A |^ + e In |e| + Ai |^ + | A-i H). 

Collect the terms of the same order of small parameter and the same 
exponents. It yields the equations for coefficients of representation (jllj] . 

dt^l +dl^l +\l\^h=Q- (12) 

dt^ ^ V +2\h\^v + VV* = {) (13) 

-S' A= /; 
-S' Ai= -dt, A -dl^ A-2\h\^A] 

S' A-,= -h'A*. 

Initial conditions for differential equations for v are obtained by matching. 
These conditions are evaluated on the resonance curve to = 0. 

v\t^=Q=u{xi,Q) + {l-i)^f{xi)] (14) 

V |i2=o =Wq (xi). (15) 

The residue part e^r{t2,xi,e) = o{e) as t2 ^ £■ This condition is deter- 
mined the domain of validity for (fTT|) . 

Formula (fT^ is connection formula for solution before and after the res- 
onance. Additional term (1 — i)y/7if{xi) leads to changing of the solution 
after passage through the resonance. 

Acknowledgments. We are grateful to I.V. Barashenkov, L.A. Kalyakin 
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